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Abstract 

Galactic magnetic fields are observed of order ~ 10~ 6 G, but their origin is 
not definitely known yet. In this paper we consider the primordial magnetic 
fields generated in the early universe and analyse their effects on the den- 
sity perturbations and the CMBR anisotropy. We assume that the random 
magnetic fields have the power law spectrum and satisfy the force-free field 
condition. The peak heights of the CMBR anisotropy are shown to be shifted 
upward depending on the magnetic field strengths relative to the no-magnetic 
field case. 
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I. INTRODUCTION 



Recently many possible generation mechanisms of primordial magnetic fields have been 
suggested to explain observed galactic magnetic fields of order ~ 10~ 6 G [[I]. The dynamo 
mechanism explains the origin of the large scale galactic magnetic field with amplification 
of a small frozen-in seed field to the observed [iG field through turbulence and differential 
rotation. The dynamo saturates when the growth enters the nonlinear regime. However the 
saturation might actually be too fast for a large scale field to form M . Without the dynamo 
mechanism, to explain the galactic fields from the primordial fields which get compressed 
when the protogalactic cloud collapses, the needed amplitude of the primordial magnetic 
fields is quite large to be of the order of 10~ 10 ~ 10~~ 9 G. Cosmological phase transitions in 
the early universe may produce magnetic seed fields. If conformal invariance is broken during 
the inflationary period, magnetic seed fields are generated And also the electroweak 
phase transition [|J || and QCD phase transition |J can generate magnetic seed fields. 
Gasperini et al. || considered generation mechanism in stringy model with broken conformal 
invariance by a dilaton field. But the field amplitudes produced by several mechanisms are 
much too weak to explain observations. 



Primordial magnetic fields may generate density perturbations || [Kj ||11|| . Tsagas and 
Barrow |T2[ |T3] considered the general relativistic density perturbations with magnetic fields. 
To treat the large scale cosmological perturbations we confront with the gauge ambiguity 
problem. It is caused by the regions larger than horizon size being causally disconnected. 



Bardeen formulated the gauge invariant method to solve the gauge ambiguity problem [|14 



The details about the gauge invariant method of cosmological perturbations can be found in 



Ref's [14 1, |15| and |[L6| . Cosmological perturbations can be classified according to how they 
transform under spatial coordinate transformations in the background spacetime; scalar, 
vector, and tensor perturbations. They relate to density, vorticity and gravitational wave 
perturbations respectively. Here in this paper we only consider scalar perturbations. From 
the observations that the magnetic field energy density is much less than the background 
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radiation energy density, we can treat it within the linear perturbation theory. 

The Big Bang Nucleosynthesis(BBN) can constrain the amplitude of magnetic fields at 
early epochs. It is argued in Ref. ]T7[ that the presence of magnetic fields affects BBN by 
changing the weak reaction rates, the electron density and the expansion rate of the universe. 
So they put constraint on the magnetic field amplitude B < 2 x 10 9 G at T = O.OlMev. 



Barrow et al. |18[ also derived an upper limit of the magnetic field amplitude at present 
Bq < 3.4 x 1CT 9 G using microwave background anisotropy created by cosmological magnetic 
fields. 

The CMB photons are polarized through the Thomson scattering of the photons and 
electrons during the decoupling time |IIJ . The upper limit on its degree of linear polarization 
large angular scales is Ap < 6 x 10~ 5 |20[. We expect that the CMBR polarization on small 
angular scales would be observed with the future experiments, MAP |33| and PLANCK [34 



If primordial magnetic fields exist at the decoupling time, they cause Faraday rotations 
which rotate the directions of polarization vectors. This effect can be imprinted on the 
cosmic background radiation and we may obtain informations on the amplitude of primordial 
magnetic fields by measuring the polarizations of the CMBR |2T|] P^[ . 

In this paper we calculate the evolution of density perturbations with the primordial 
magnetic fields which have power law spectrum. We do not concern ourselves with the 
details of generation mechanism of magnetic seed fields, but assume that sometime during 
radiation dominated era, large scale magnetic fields are generated instantaneously. We then 
investigate how they affect the temperature anisotropy and polarization of the CMBR using 
various spectral indices and field strengths of the magnetic field. 

In section [D], we derive, using the gauge invariant variable, the density perturbation 
equations with magnetic fields present. The equations are solved numerically and the effect 
of magnetic fields on the temperature anisotropy and polarizations of the CMBR are shown 
in section III. Finally we summarize the results in section IV. 
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II. DENSITY PERTURBATIONS WITH MAGNETIC FIELDS 



In this section we consider the background space is homogeneous and isotropic. Cosmo- 
logical perturbations are classified as scalar, vector and tensor perturbations according to 
how they transform under spatial coordinate transformations in the background spacetime. 
We will treat here only scalar perturbations which are related to density perturbations. 
In the longitudinal gauge (conformal Newtonian gauge) the metric, including the scalar 
perturbations, is written by [ITH] , 



ds 2 = a(r]Y [-(1 + 2^)d7] 2 + (1 + 2^ ij dx l dx 3 ) (1) 
where rj is the conformal time defined by dt = a(rj)dr] and 7^ is the spatial metric tensor. 



^ and $ are related to the gauge invariant quantities $^ and of Bardeen fllj] and the 



gauge invariant potentials \1/ and $ of Kodama & Sasaki fL5 |. The physical meaning of $ 
and \l/ are the curvature perturbation and Newtonian gravitational potential respectively. 
The Maxwell equations have the form 

F" v . u = J» (2) 

Fp,v\p Fup;)i ~l~ Fpfj,;u 0, (3) 

where F pv is the second-order antisymmetric Maxwell tensor which represent the electro- 
magnetic field and J M is the 4- vector current which generates the electromagnetic field. The 



Maxwell tensor splits into the electric and magnetic 4- vector [|12[], defined by 



E» = F^u", (4) 
B, = ^, upX u»Fr\ (5) 

where e pup \ is Levi-Civita tensor and is the fluid four-velocity. The background value of 
is taken to be = (1/a, 0,0,0). Then the electric and magnetic 4-vectors are purely 
spatial, i.e. E^u p = and B^u p = 0, so we denote the spatial components E % and B % by E 
and B. 



The generalized covariant Ohm's law is 

J» + J v u v u^ = oF^u v (6) 

where a represents conductivity of the medium. The spatial components of Eq. (|6|) are 
reduced to J = crE where J is the spatial component of J M . Assuming infinite conductivity 
of the medium in the Universe ||, we neglect the electric field so that E = 0. 
Now we can reduce Eq. (Q), using the magnetic field 3- vector B(??, x), to 

«g?U V-B^O, (7) 

where V is the covariant differentiation with respect to 7^. In this work, we consider only 
the case where — $ij f° r simplicity. From the first of Eq. (J7J), we get B(?7,x) oc a 3 . 
The magnetic field energy density, \B 2 {= ^B^B^ = \a L Yh=i B l B l ), evolves the same as 
the radiation energy density ~ or 4 . The dimensionless quantity r is introduced in Ref. f3[ 
defined by r = -0- which is the ratio of magnetic field energy density to the background 
radiation energy density. It is approximately constant at all early history of the Universe. 
Total energy momentum tensor is decomposed by 

rpflU rp(fluid)flV _|_ rp[em)fJ,U ^g^ 

where the electromagnetic energy momentum tensor and fluid energy momentum tensor are 
given by 

T (em)nu = p^puX _ \g^ Fp(jF P^ (9) 

and 

T {fiuid)„u = ( p + py» u v + Pg ^u ( 1Q ) 

We here treat the matter as perfect fluid to neglect the anisotropic pressure perturbations 
and consider only adiabatic perturbations to neglect the entropy perturbations. The lin- 
earized perturbation equations are obtained from the Einstein equations up to first order, 
and are written as follows; 



3 q, - 3-V + (V 2 - 3K)$ = -AuGa 2 p ( A + -—) , 




(11) 



a' 

Vi(-tf - $') = -4vrGa 2 (p + P)^ 



(12) 



a 




(14) 



(13) 



where the prime denotes derivative with respect to the conformal time 77, and lT em ^j = 



gauge invariant density contrast and velocity perturbation and c s is the sound velocity V 2 
is the Laplace-Beltrami operator whose eigenvalue is —A; 2 . Eq. (|H|) is the trace part of the 
spatial component of the perturbed energy-momentum tensor and Eq. ( |i~4"|) is the traceless 
part. 

To write down the perturbation equations for a given wave mode, k, we define the 
Fourier transform of the perturbed quantities and random magnetic fields. In this paper we 
consider density perturbations in flat space, K = 0, so the spatial dependence of the Fourier 
transform is just the plane wave, e* k ' x (= e lkiX *), 



and also v(x), $(x),\l/(x) and B(x) are defined similarly. We assume that the force-free 
field condition (B x V x B = 0) is satisfied to treat magnetic field. Then Eq.'s (|TTD ~ flUD 
can be written by 



\$ l jB 2 — B l Bj corresponds to the magnetic field anisotropic pressure. A and v $ are the 




A(k) 



(15) 




(16) 



a' 

fe(— * - <£>') = -AixGa 2 {p + P)v 



(17) 




(18) 
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fc 2 ($ + tf) = -87rGa- 2 Q^W) ( 19 ) 
where F(k) is defined by 

F(k) = J d 3 k'B l (k')B l (k - k'). (20) 

which represents the spectral dependence of magnetic fields. The fact that the magnetic 
field energy density decays as ~ a -4 is used. In the Appendix A, we calculate the Fourier 
transform of the magnetic field anisotropic pressure using force-free field conditions and 
derived the expression F(k),Eq. (^0[). 

To investigate the spectral dependence of perturbed quantities, we need to take ensemble 
average of \F(k)\ 2 due to random magnetic field. For a homogeneous and isotropic random 



magnetic field, B(k) satisfy the relation Jl0| [23] 



< B\k)B^) >= 5 3 (k - k') - ^ (21) 

and then 

< Bl >= J d 3 kB 2 {k) (22) 

where angular brackets denote a statistical average over an ensemble of possible magnetic 
field configurations and < B% > x l 2 is the average field strength observed today. 
From Eq. (Eg) 



|F(A;)| 2 = J d 3 k' d 3 k" B(k') ■ B(k - k')B(k")* ■ B(k - k")* (23) 

Taking ensemble average of the both sides, we obtain 

< |F(A;)| 2 > = J d 3 k' d 3 k" < B(k') • B(k - k')B(k")* • B(k - k")* > 

= J d 3 k' d 3 k" [< B\k - k')5* m (k - k") >< S z (k')K(k") > 

+ < B\k - k')B* m {k") >< ^(k')5* m (k - k") >] (24) 



Using Eq. (p2l|) , and integrating the delta function, we can get 



. ( {(k - k') • k'} 2 ] B 2 (\k - k'|) B 2 (k') 



|k-k'| 2 A;' 2 



< |F(A;)| 2 > = 2 J d 3 k' 

/ dk'k' q+2 d^[k 2 (l + fx 2 ) + 2k' 2 -Akk'ti] 
Jo J-i 



?7T 2 



x[k 2 + k' 2 -2kk'i2]i- 1 (25) 

where k • k' = kk'fi, \x = cos# and S 3 (k = 0) = V/(27i) 3 is used. V is the volume factor. We 
used the power-law spectrum 

B 2 {k) = Ak q . (26) 
We define the average field on scale A by 

< B 2 > x = [ d 3 kB 2 (k)exp ( . (27) 



Then we can determine the coefficient A in Eq. 

1 A«+ 3 



A= -^2^mr(n+l) <B2 >X - (28) 



To determine k max we use the argument in Ref. |TT]. Small scales reach nonlinear variance 
(A > 1) earlier than large scales, and the first scales to become nonlinear have k ~ 2k max . If 
we choose the scale that corresponds to the formation of galaxies at z n i = 6, k max = irMpc -1 . 
If instead we require that magnetic fields form clusters of galaxies at z n i = 1, that would 
correspond to k max = -Mpc' 1 . As we shall see below, the CMBR anisotropy much depends 
on the k max value. Treating k/k max as a small parameter, the leading term of the result of 
integrating Eq. (|25|) is 

< |F(A;)| 2 >~ ^X-A 2 ^—k 2q t!. (29) 
1 V n (2tt) 3 2g + 3 max V ' 

For the case of q < —3, the integration of Eq. (|25"[) diverge as k' — > 0. So we only consider 

q > -3. 



III. CMBR ANISOTROPY 



Cosmological density perturbations cause the temperature fluctuations when the pho- 
tons decouple from the thermal bath at last scattering surface. Furthermore small metric 
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perturbations induce bulk velocities of the fluid, and the resulting anisotropics in the photon 
distribution will induce polarization when the photons scatter off charged particles (Thom- 
son scattering) |2T|] . After decoupling, the photons freely propagate along the geodesies, and 
any polarization produced through the Thomson scattering will remain fixed. The evolution 
of the CMB anisotropics is described by a set of radiative transfer equations. Temperature 
and polarization anisotropies are expressed in terms of Stokes parameters I,Q,U and V. 
The parameter I gives the radiation intensity which is positive definite, Q and U represent 
the linearly polarized light and V describes the circular polarization. The degree of linear 
polarization Ap is defined in terms of Q and U, Ap = (Q 2 + U 2 ) 1 ^ 2 , and the temperature 
anisotropy is denoted by A T (= AI). A T and Ap can be expanded in multipole moments 
defined such that A (77, k,fi) = J2i(2l + l)Az(r/, k)Pi(fi) where Pi is the Legendre polynomial 
of order I and \x is the cosine angle between the wave vector and the direction of observation. 
Their evolution equations are given by [p7j] 

A T + ik^{A T + #) = -$'_ K '[A T - A To - fj,v b + \P2{p)S p ], (30) 
A' P + ik/iAp = -k'[A p Pz(ti))S P ] (31) 

where Sp = — At 2 — Ap 2 +Ap . k' is the differential optical depth defined by k' = x e n e aTa/ clq 
with x e the ionization fraction, n e the electron number density and <jt the Thomson scatter- 
ing cross section. The Thomson scattering cannot produce any net circular polarization ITS 



and thus we expect V — for the microwave background. v b is the velocity perturbations 
of the baryon component which is affected by the existence of primordial random magnetic 
fields. Metric perturbations $ and which evolve according to the equations in section |TI 
under the influence of random magnetic fields, act as the source terms in Eq. ( |30"|) which 
governs the evolution of the temperature anisotropy. 

Eq.'s (|30D and ( j31~l) are formally integrated to yield p7 



ArM = r dve ix »g(v)[^T (v) + Mv) ~ ^(^Spiv)} + T drje^e-^^' - $'), (32) 
Jo 2 Jo 

A P ( m ) = jT drje^gW-il - P 2 ^)}Sp(r,), (33) 



where 



g(r]) = K'e~ K ^ (34) 



is the visibility function and 



K (Vo, V) = K dr] (35) 
J v 

is the optical depth to photons emitted at the conformal time rj. The visibility function 
represents the probability that a photon observed at 770 last scattered within dr] of a given 
rj. For the standard recombination this function has a sharp peak at the conformal time of 
decoupling i] D . And x = k(r] — rf). 

Under a clockwise rotation in the plane perpendicular to the direction of observation, n, 
the temperature is invariant while Q and U transform as 

Q' = gcos2^ + [/sin2V>, 
U' = -g S in2^ + Ucos2ip, (36) 

or 

(Q±iU)' = e^(Q±iU) (37) 
where ip is the rotation angle. Therefore the quantities can be expanded in terms of the 



spin-2 spherical harmonics [28 



(Q±iU)(h) = J2 a ±2M±2Yim(n) (38) 

im 

where ±2Y l m (n) is the spin-2 spherical harmonics whose properites are summarzied briefly 
in Appendix B. The expansion coefficients are 

a±2 M = J dQ± 2 Y l * m (Q ± iU)(n). (39) 



In Ref. [PS| |, the authors introduce the following linear combinations of a±2,/m to circumvent 
the difficulty that the Stokes parameter are not invariant under rotations; 
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aE,lm — —-^{ a 2,lm + < l ~2,lm) 
1 

0-B,lm = — (d2,lm — a -2,lm)- (40) 

These newly denned variables are expanded in terms of ordinary spherical harmonics, Yi m , 

E(n) = ^a E ,im Y irn(n), 

I in 

B(h) = ^a Biim F im (n). (41) 
i in 

The spin-zero spherical harmonics, Y/ m , is free from the ambiguity with the rotation of 
the coordinate system, and therefore E and B are rotationally invariant quantities. The E- 
mode has (—1)' parity and the B mode (— parity in analogy with electric and magnetic 
fields. Scalar perturbations generate only the E mode of the polarizations []3D|]. The power 
spectra of temperature and polarization anisotropies are defined as Cn =< |«r,Zm| 2 > for 
At = Y,irn a T,imYim and analogously for C E i- So if we get the evolution of the temperature 
and polarization anisotropy amplitude from Eq. ([32]) and (|33|) , the amplitudes for each mode 
of power spectra are given by 

C T j = (4vr) 2 j k 2 dkP s (k)[A T ,t(k)] 2 , (42) 
C E)l = (4tt) 2 J k 2 dkP s (k) [A Ejl (k)} 2 , (43) 
Cci = (47r) 2 | k 2 dkP s (k)A Tl (k)A El (k) (44) 



where Ps{k) is the initial power spectrum and Ati and A E i are given by p8| , 
rvo 

A Tl (k)= / dr}S T (k,ri)ji(x), (45) 
J 

(1 + 2Y n° 

A Ei(k) = , I dtiSsfcTfijtx), (46) 



\ (I -2)1 Jo 

S T (k, „) = ,(A ro + * - g - * - f ) + .-<* + *) - 9 >q + - M (47) 
fc(M) = |£. (48) 
We here concern ourselves with the flat CDM universe with adiabatic initial conditions. 



We use the CMBFAST code |26| to calculate numerically the CMBR anisotropy. During this 
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calculations we put h (Hubble constant divided by lOOkm/sec/Mpc) = 0.5 and assume 3 
species of massless neutrinos. In Fig. 1, we plot the angular power spectrum of temperature 
fluctuations 1(1+1)Cti with the magnetic field strengths, 3x 10~ 8 G, 5x 10~ 8 G and 7x 10~ 8 G, 
for a given magnetic field spectrum index, q — 1. Observed amplitude of galactic magnetic 
fields is order of ~ 10~ 6 G. The BBN can constrain the amplitude of magnetic fields, B < 



10 G jl7| , and also derived an upper limit of the magnetic field amplitude B < 10 G using 



the CMB anisotropy |18||. Another constraint on magnetic field intensity can be obtained 
from r (= B 2 /2p r ) < Ah where Ah is the horizon crossing amplitude. COBE 4-yr data 
gives Ah ~ 10~ 5 on large angular scales. To consider magnetic fields of order of ~ 10~ 8 G 
does not violate too much current bounds on magnetic field amplitudes by the observational 
and theoretical considerations. The figure shows that the spectral curves of the CMBR 
temperature anisotropy are shifted upward with increasing magnetic field strengths. We 
can conclude from the numerical calculations that the presence of the magnetic fields which 
have field strength of order < 10 _9 G at present doesn't significantly affect the temperature 
fluctuations. The density perturbations with magnetic fields of order of 10 _8 G result in the 
deviations of the angular power spectrum, Cti, of up to 14%. 

The i?-polarization spectrum, 1(1 + 1)Cei, and temperature-polarization correlation spec- 
trum, 1(1 + l)Cci, are shown in Fig.'s 2 and 3 for q = 1. Also in these figures we can see 
that the spectrum curves are shifted upward with increasing magnetic field strengths rela- 
tive to the non-magnetic case. The current bound on the degrees of linear polarizations of 
the CMBR on large angular scales is A P < 6 x 10~ 5 [pOj] . As we discussed in the previous 
section, we plot the temperature anisotropy with k max = n/2Mpc~ l and nMpc~ l in Fig. 4. 
In this figure, we can see that there is a strong dependence of the spectrum curves on the 
cutoff k max . 

In Fig. 5 we plot the temperature anisotropy with the spectral index of magnetic field 
q — 1,2 and 3 for B\ = 5 x 10~ 8 G with A = O.l/i^Mpc. The spectrum curves are nearly 
independent of the spectral index. We probe the vicinity of the acoustic oscillation peak, 
/ ~ 200, to investigate the dependence of spectral index more closely. The result is that 
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the spectrum curves are shifted downward with the increasing spectral index. In Ref. 
p5| recently, the authors derive an expression for the angular power spectrum of CMBR 
anisotropies due to gravity waves generated by a stochastic magnetic field. They show 
that, for n > —3/2 (n is magnetic field spectral index in their notations), the induced C\ 
spectrum from gravity waves is independent of n, but only the amplitude depends on the 
spectral index, l 2 Ci ~ (A/c ma:r ) 2n+3 / 3 . They also derive an upper bound of B\ for n > —3/2 
and A = O.l/i^Mpc 

B x < 9.5 x l(TV a37ri G. (49) 

Here we don't consider the Faraday rotation due to the magnetic field which can change 
the polarization spectrum because we restrict our calculations in linear perturbation theory. 



The authors in Ref's ||21|| , p2| and |32| studied the effect on the CMBR anisotropy with the 
uniform primordial magnetic field causing Faraday rotations in the homogeneous background 
universe. They argued that the presence of magnetic fields depolarize the CMBR anisotropy 



2"2|] and proposed that the temperature and B mode polarization correlation which are 



generated by Faraday rotations can constrain the magnetic field 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper we consider the density perturbations with primordial magnetic fields 
present using gauge invariant formalism. While magnetic field generation mechanism is 
not yet known, we assume that the magnetic fields smear out all over the universe randomly 



in the radiation dominated era. Using the CMBFAST code |26| we solve numerically the cou- 
pled density perturbation equations for flat CDM universe with adiabatic initial conditions. 
We investigate the CMBR anisotropies for the magnetic field permeated universe. With the 
temperature anisotropy spectrum we cannot fully determine the cosmological parameters 
and the informations about the cosmological perturbations. CMB photons are polarized 
due to Thomson scattering during the decoupling time. Small cosmological perturbations 
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induce the polarization at last scattering surface through the Thomson scattering. The lin- 
ear polarization relates to the quadrupole anisotropy in the photons. So if we investigate 
the polarization as well as the temperature anisotropy, we can get enough information at 
the time of decoupling. The information can constrain the cosmological parameters and the 
cosmological perturbations. Next we study the effect of the random magnetic field on the 
CMBR temperature anisotropics and polarizations. We consider the several scale magnetic 
fields with the assumption of power law magnetic spectrum. To get the polarization power 
spectra, we use the rotation invariant scalar quantities E and B which are introduced in 



Ref. [^8j. B vanishes for scalar perturbations. 

For a given spectral index the temperature anisotropy and polarization spectrum are 
shifted upward with increasing magnetic field strengths. The density perturbations with 
magnetic fields of order of 10~ 8 G result in the deviations of temperature anisotropy power 
spectrum of up to 14%. The fluctuations of this order due to the primordial magnetic field 
are sufficiently large to be observed in future satellite experiments. Further, if B\ <~ 10 _9 G, 
magnetic field energy density does not affect noticeably the CMBR anisotropy. We assume 
that the magnetic fields have the power-law spectrum. The spectrum curves are nearly 
independent of the magnetic field spectral index. 

Here we assume that the magnetic field energy density evolves as ~ a -4 . But in the 
early era, when the magnetic fields are generated, their evolution behaviors may be different 
depending on generation mechanism. If so, temperature fluctuations due to magnetic field 
may be shown. 



In the early next century, the new satellite experiments, MAP [B3] and PLANCK [34 



will be set forth with better accuracy than COBE satellite. They are expected to detect the 
imprint of the polarization as well as gravitational wave. If it is possible, we can constrain 
the magnetic field strength and the spectral index and get the hint about the magnetic field 
generation mechanism. 
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APPENDIX A: CALCULATION OF TRACELESS PART OF PRESSURE 

PERTURBATION 

Here we calculate the Fourier transform of random magnetic fields with the assumption 
of for the force-free magnetic field (B x V x B = 0) condition. From the perturbed Einstein 
equations, we can write the traceless part of pressure perturbations as follows in the real 

space, 

(V l Vj - + = -87rGa 2 n( em ^, (Al) 

3 

where Yl^ em ) 1 ■ = |5*5 2 — B l Bj. Using the Fourier transform of B(x), n <em ) l j(x) can be 
written by 

n (em)^ x ) = J ( fkd z k'[U)B l i}i)B l {\i - k') - 5 J (k)^(k - k')]e lk ' x (A2) 

where we omit the time dependence for brevity. Then we differentiate n^ em ^(x) with respect 
to x % to get, 

V;n( em ^-(x) = J d! i k<fk'ik i [-5)B l {\i)B l (\i - k') - 5 l (k)^(k - k')]e ik - x (A3) 
We can assume 

J ^k'B'i^Bjik - k') = {A5) + B^)k 2 F(k) (A4) 

which is obvious from the fact that tensorial component of scalar perturbations is split into 
the trace and traceless part. 

Next, differentiating B l Bj with respect x % yields 

V i (S i S j ) = (B x V x B)j + &VjBi. (A5) 
15 



The first part of the right-hand side is the magnetic force due to the current density(J = 
V x B), and we neglect this term assuming that the force-free condition is satisfied in the 
early universe. We take Fourier transform of the lT em ) 4 ,,- and then again differentiate with 
respect to x l using the force-free field condition to obtain, 

ViII (em ^(x) = Jd 3 kJ ^k'ik^B^k^B^k - k') - ^(k')^(k - k')]e ik ' x (A6) 

Comparing Eq. (A3) and Eq. (IMS), we can find the relations 

A + B = - 
2 

F(k) = J d 3 k'B l ik!)Biik - k') (A7) 

Then Eq. (|A1|) is written by 

[ d 3 k(-5) - ^)k 2 ($(k) + ^(k))e ik ' x = -8nGa 2 [ d 3 k(-S) - A5) - B^)F{k)e i ^ (A8 
J 3 J fir J 3 J J k 

Further we only treat the scalar density perturbations and traceless component of pressure 
perturbations, so the right hand side should be proportional to (|5j — ^l 1 ). Then we find 
that A must have value of \. Finally, we can write the traceless part for a given mode k, 

A; 2 ($ + ^) = -8nGa 2 (^F(k)) (A9) 



APPENDIX B: SPIN-S HARMONICS 



In this appendix we summarize the definition of spin-s function and the property of spin-s 



harmonics. We mainly refer to Ref.'s p8fl and [29 



A function s f(9, 0) defined on the sphere is said to have spin s if under a right-handed 



rotation of (ei,e*2) by an angle ip it transforms as s f'(9,4 
tions can be expanded in spin-s spherical harmonics, s Yi m (t 
orthonormal basis. The spin-s harmonics are expressed as 



■ e s ls ^f(9,(p). A spin-s func- 
), which form a complete and 



(l + m)\(l-m)\2l + l 
(l + s)\(l-s)\ 



1/2 



sm 2/ (^/2)E 



(l 




' l + s 


\ 


\ 


r ) 


^r + s — 


m ^ 



x (-l) / - r '- s+m cot 2r+s - m (^/2). 



(Bl) 
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These set of functions satisfy the conjugatation, completeness and orthogonality relations: 



2?r 



dcos9 s Y i r m ,(9,(p) s Y lm (9,(f)) = 5 V1 5 V 



J2 Y; m (9, 4>) s Y lm {9\ 0') = 5(0 - 0')<5(cos# - coaff). 

lm s 

Finally the harmonics are related to the ordinary spherical harmonics as 
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FIGURES 
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FIG. 1. The angular power spectrum of temperature fluctuations with the magnetic field 

strength B A = 3x 10 _8 G, 5 x 1(T 8 G, 7 x 10 _8 G with spectral index q = 1 for A = 0.1/t _1 Mpc. 
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FIG. 2. The E mode polarization spectrum /(/ + 1)Cei for q = 1. 
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FIG. 3. Temperature and polarization cross-correlations, 1(1 + l)C c i, for q = 1. 
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FIG. 4. The angular power spectrum of temperature fluctuations with the value of the magni- 
tude of cut-off wavevector k max = nMpc" 1 , n/2Mpc~ l for B\ = 5 x 10~ 8 G for A = 0.1/t _1 Mpc. 
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